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We report on the reinforcement of superconductivity in a system consisting of a narrow supercon- 
ducting wire weakly coupled to a diffusive metallic Him. We analyze the effective phase-only action 
of the system by a perturbative renormalization-group and a self-consistent variational approach to 
obtain the critical points and phases at T = 0. We predict a quantum phase transition towards a 
superconducting phase with long-range order as a function of the wire stiffness and coupling to the 
metal. We discuss implications for the DC resistivity of the wire. 
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I. INTRODUCTION 

The interplay between fluctuation and dissipation phe- 
nomena in quantum systems is presently under intensive 
research. Fluctuations are particularly strong in low di- 
mensions, as reflected by the lack of long-range order 
in ID-systems with short-range interactions^!. On the 
other hand, dissipation counteracts fluctuations effects, 
decreasing the lower critical dimensiorP^. 

Some physical realizations of dissipative low- 
dimensional systems are the well-known resistively 
shunted Josephson junctions arrays, where the effect of 
local ohmic dissipation has been intensively studiecP^, 
superconducting grains embedded in metallic films 12 " 15 
and Luttinger liquids coupled to dissipative bathiP^HUl 

Narrow superconducting wires with diameter d <C Co 
(where £o is the bulk superconducting coherence length) 
are low-dimensional systems in which strong fluctuations 
of the order parameter affect low-temperature properties. 

It was originally suggested by Little^, and subse- 
quently discussed by Langer and Ambegaokar (LA)G^ 
and McCumber and Halperin (MH) 2 3 , that resistivity 
q (T) in thin wires would be finite for all temperatures 
below the bulk critical temperature T c . Thermal fluc- 
tuations cause the magnitude of the order parameter to 
temporarily vanish at some point along the wire, allowing 
its phase to slip by 2ir (the so-called thermally activated 
phase-slips) and dissipate through the Josephson relation 
V = h/2e AO, where A8 is the phase difference across the 
wire. 

According to the LA-MH theory, thermal fluctuations 
induce a resistivity g(T) ~ f2 (T) exp [— AFq/T], where 
AFo is the Ginzburg-Landau free-energy barrier between 
different current-carrying states in the wire, and Q (T) is 
an algebraically decreasing function of T. However, de- 
viations from the LA-MH theor y in t he regime T«T C 
were first observed by Giordano 1 22 * 23 ! and more recently 
by other experimental group d 24 " 28 !, leading to the conclu- 
sion that for very thin wires at low temperatures current- 
decay was produced by macroscopic quantum tunneling 



of the phase of the order parameter through the same 
free-energy barriers (the so-called quantum phase-slips), 
leading to a much weaker dependence of the resistivity 
on T. 

Moreover, it is believed that the destruction of the su- 
perconducting state in very thin wires occurs through 
the proli feration o f quantum phase slips/anti phase 
slips pair d 23 * 24 * 27 * 29 * ^, in what constitutes the quantum 
analog in 1+1 dimensions to the classical Berezinskii- 
Kosterlitz-Thoulesd^ transition in two space dimensions. 




Figure 1: (color online) Representation of the system. At T <C 
T c one-particle hopping is suppressed by the BCS gap-energy 
Ao. At next order in the hopping process, Cooper pairs can 
tunnel into the metal and propagate coherently in the metal, 
generating an effective coupling ~ cos [6 (r) — 6 (r')) in the 
wire. 

Contrary to other ID systems such as dissipative ohmic 
Josephson junction arrays, isolated thin wires d o not 
present significant sources of dissipation at T -C rj 30 l 33 [ 
However, additional sources might be provided by a cou- 
pling to the environment, a possibility which has hardly 
been explored yet. Although general theoretical frame- 
works have been pro posed to describe superconductor- 
normal (SN) junction d 13 * 35 * 36 !, recent advances in super- 
conducting nanowires fabrication techniques call for a 
more detailed analysi s of the effects of coupling to general 
dissipation source9^ h 24 | 27 | 31 l 

In this article we focus on the effect of weakly coupling 
a superconducting wire to a diffusive 2D normal metal. 
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We show how the induced dissipation stabilizes supercon- 
ductive long-range order at T = despite the ID nature 
of the wire. At finite T, the effect of dissipation are man- 
ifested in an increase of the superconductive stiffness of 
the wire. 

The article is organized as follows: in Section [TT] we 
derive the effective low-energy phase-only action of the 
coupled system. Section [TTT| is devoted to the analysis of 
this model within a pertubative renormalization group 
and a self-consistent harmonic approximation and discuss 
implications for the DC resistivity. Finally in Section |TV| 
we discuss the main physical consequences of our results 
and summarize them. 



II. MODEL 

We analyze the system depicted in Fig. [T] which rep- 
resents a clean superconducting wire of length L and 
lateral dimensions d <^ ^q, weakly coupled to a diffu- 
sive 2D metal. In the following we use the convention 
h = kg = 1 . We begin our description with the action of 
the microscopic BCS Hamiltonian for the isolated wire 



Sy, 



J\t J d 3 R (3t-^ + - 



vy v] [V^ g 
2m 



+ U J dr I d*R VtV-^f 



(1) 



Here the fermionic field tp a = tp^ (R, r) describes an 
electron in the wire with spin projection a at position 
R = (x, y, z) and imaginary-time r. The chemical poten- 
tial /j, = hp /2m is the Fermi energy in the normal state 
and the local attractive interaction U < is responsible 
for pairing at T < T c . 

Assuming for simplicity that the coupling to the metal- 
lic film takes place along the line (x, 0, 0) in the wire, the 
coupling term is described by 



S± 



tj_ I dr 
h 



dx 



[Vv (x, r) y] a (x, r) + h.c] , (2) 



where the fermionic field r) a (r,T) represents an elec- 
tron at position r = (x,y) in the film. Here the com- 
pact notation ip a (x,T) = ip a (R,r)\ y=z=0 , r] a (x,T) = 
r) a (r,T)\ y _ has been used. While certainly more re- 
alistic models for the coupling, which take into ac- 
count geometr ical details of the SN junction have been 
studied^- 3 * 35 * 36 !, the main physics which is of interest to us 
is already captured by Eq. ^ . 

Electronic motion in the metallic film is described by 
the non-interacting action 



S 2 d = [ dr f d 2 r S~] {f]o (d T - ^2d) 
Jo J 



2m 



+ V (r) fj a (r) r, a (r) \ , (3) 



where V (r) is the (static) disorder potential. We assume 
weak enough disorder, such that the localization length 
in the film is £i oc 3> L. allowing us to neglect strong 
localization effects. 

For one given realization of the disorder potential 
V (r), the effective action in the wire is obtained by inte- 
grating the electronic degrees of freedom in the metallic 
film 



^diss 

p 



^diss = — 1± J drdr J dxdx x 

(7 

where gm (r,r';r — t j ) is the Green's function in the 
film. Note that the spin index has been dropped using the 
SU (2) symmetry of the problem and that we used the 
notation g 2D (x,x';r- r') = g 2D (r,r';r- t')\ 

Since the disorder potential breaks the original trans- 
lation invariance in the wire, an average over different 
configurations of the disorder is needed to restore it. Let 
us define the partition function of the systems for one 
disorder realization as 



Z[V} = Jv fofl e 



Sw-Sd 



Assuming for convenience that V (r) is gaussian- 
distributed 

Sd= 2^ j* T V2(rh 

we can formally perform the average over different disor- 
der configurations as 



Z 



J V[V] e~ Sd Z[V} 
JV[V] e~ s * ' 



Expansion of Z [V] in powers of t± allows us to obtain 
an explicit form of the partition function Z 



J V [V] e' 8 * 

r 00 i 

«* n 



(4) 



n=0 



The low-energy effective action of this model is 
obtained by introducing Hubbard-Stratonovich fields 
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A (R, t) , A* (R, t) in the pa rticle-part icle channel to de- 
couple the quartic term in ,g w l 3 l 33 l 3 Z1281 After integration 
of the fermionic degrees of freedom in the wire, the action 
reads 



where 



AT) 



S w [A, A] = -Tr Ing- 1 - i J dr J d 3 R |A (R, r)\ 2 



(5) 



where 

'd T -n-^ A(R,r) 

A(R,r) a T + / i+^ 

and where the Nambu notation 

_ ff t (R,r) 



*(R,t) 



^i{R,r) 



is implicit. The Green's function in the wire formally 
reads 



9w (R, t) 



g(R,r) f (R, r) 
J{R,t) g{R,r) 



where g (R, r) = (T T i >{R,r)ip (0)) and g(R,r) = 
(T T ijj (R, t) ip (0)) denote respectively the particle and 
hole propagators in the wire, while / (R, r) = 
(T t ^(H,t)^(0)), f(R,r) = (T t $(R,t)$(0)) are the 
anomalous ones^l. 

For very narrow wires with diameter d <C £o at low 
energies, the dependence of the fields A (R, r) , A (i?, t) 
on transverse dimensions can be neglected, reducing to 
A (R, t) — > A (x) where the compact notation x = (x, r) 
has been used. Moreover, at T <§; Tmf (where Tmf is the 
mean- field critical temperature of Eq. (|5|) and neglecting 
amplitude fluctuations, the dynamical state of the wire 
is characterized by A (a;) = A Q e ie ( x \ where the quan- 
tity Ao corresponds to the (temperature-dependent) BCS 
energy-gap and 9 (x) is the space- and time-dependent 
phase of the macroscopic BCS wavefunction. 

The derivation of the phase-only action in the isolated 
wire (i.e., the first term in the expansion of Eq. Q) 
is obtained by the means of an expansion in Gaussian 
fluctuations in the gradients of the field 9 {x) around the 
BCS sadd le-point, an d takes the form of a Luttinger liq- 
uid actioiPasmiaini 



So 



dx 



(6) 



Here II = II (x) is the momentum canonically con- 
jugate to 6(x). formally defined through the relation 
[8 (x) , II (x 1 )] — iS (x — x') and representing fluctuations 
in the density of Cooper pairs at point x. The operator 
V denotes derivation with respect to the spatial coordi- 
nate x. The Luttinger liquid parameters u and K are 
defined atP^ 



!A w n 8 (T) 
Amn (T) 



K = 2tt 



A w n s (T) K (T) 
4m 



, is the superconducting stiffness of the wire 
(with n s (T) the 3D density of electrons in the condensate 
and m their mass) , A w is the cross-sectional area of the 
wire and k (T) is the compressibility (cf . Ref. [30] for 
more details), u corresponds to the velocity of the plasma 
(Mooij-Schi^irP2) mode. In the following, we assume the 
wire to be in the thermodynamic limit L 3> Lt = u/T. 

The next terms in the expansion of the partition func- 
tion (Eq. Q) provide the effects of the coupling to 
the metallic film. At order O (t 2 ^ and at low temper- 
atures (T <C Tmf), the transfer of individual electrons is 
strongly forbidden by the energy gap Ao, giving a prob- 
ability ~ e~ A °/ T for such a charge transfer channel. 

The most relevant contribution of the coupling to 
the metallic film appears at order 0(t\), and cor- 
responds totheAndreev reflection occurring at SN 
interface d 12 * 13 * 35 * 36 * 42 !. This contribution physically repre- 
sents processes in which paired electrons (for which there 
is no energy cost) are effectively transferred from the wire 
to the film and viceversaP^l 



S A 



1 



(•Sdiss) 



= 2*1 



J dx" f{x") 



J dxdx'P c (x - x 1 ) cos [9 (x) - 9 (x 1 )] , 



(7) 



where f (x) = f (R, <r)| y=z=0 . The kernel P c (x) = 
P c (r, t)L_ is the cooperon propagator in the diffusive 
film, defined afPHU 

P c (r -r',T- t') = (g 2D (r, r', r - r') g 2 u (r, r', r - r')) d , 

(8) 

representing the probability to find a coherent electron 
pair traveling a distance |r — r'\ in the interval t — r' 
through the disordered filrrPH (see Fig. [lj. The dif- 
fusive propagation of this electron pair remains phase- 
coherent over a length £y (assumed 3> ^o) which de- 
pends on T, magnetic field and the strength of Coulomb 
interactions 43 ' 1 . In the absence of the latter, £jv (T) — 

y^P, which leads to important non-local coupling effects 
at low enough temperatures. 

Explicit evaluation of Eq. Q for a diffusive 2D 
metal, assuming a Fermi-liquid description, yields (see 
Appendix |Al 



(9) 



where we have defined 

exp ( - j- 
P c (r,r)=Re> ' ^ 



4Df I ( ir 

rfo, — 



Here T (a, z) is the incomplete gamma function and f = 
t + ir e , with r e the elastic lifetime of electrons in the 
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diffusive filrrP. Eq. ^ is a valid expression for r ^> r e 
and x 3> lei where l e is the elastic mean-free path in the 
film. In what follows, we set y — and consider the 
kernel as depending only on the coordinate x. 

The coherence length £jv (T) separates two distance 
regimes of interest: 

(a) the local regime x 3> <t/v (T), where the cooperon 
can be considered local in space, reducing to 



(*)«^Mh/ J/)r 



consistent with the expression for the Andreev conduc- 
tance in Refs. [35] Introducing the notation q = (k,LJ m ), 
where uj m is the bosonic Matsubara frequency uj m = 
2irmT, the approximated Fourier transform (neglecting 
the logarithm) is independent of k for k < ^ and reads 



Pc (q) « 2£ 



N 



1 



in the limit q — >• 0. 

f&J TTie non-local regime of distances a; < £p? (T), 
where Eq. Q can be approximated as 



P c (x) W (4D) a 
x 4 + (4Dr) 2 



(10) 



Luttinger parameters if and u of Eq. ([6| (e.g., diffu- 
son propagatorS2l) . Although these contributions do not 
change the physics at a qualitative level, their effect is 
relevant for the comparison with real systems. A micro- 
scopic study of the dependence of K and u on the hop- 
ping tj_) as we U as further renormalization arising from 
Coulomb interactions between the wire and the film, is 
beyond the scope of the present article and will be given 
elsewhere^. In the following we assume that the Lut- 
tinger parameters appearing in Eq. (|6| already include 
all these corrections. Note also that the coupling to the 
metal modifies the bare value of Ao through the well- 
known proximity-effect, in which the diffusion of normal 
electrons in the superconductor produce a lowering of 
T c E2l. However, since this is a small effect of order O (tj_) 
and in addition we assume r « T c . this effect is irrele- 
vant to our description and can be effectively taken into 
account in renormalized values of Ao and T c . 

So far we have not included the effects of topological 
defects (phase-slips) in the wire. As discussed in Sec. [T] 
these topological excitations produce finite resistivity at 
T<T C and are believed to be the origi n of destruction of 
superconductivity in narro w wir e d 30 * 3 ^ and in dissipative 
Josepshon junctions arraya^EH it can be showrP that 
defining a field <j>{x), such that W(j)(x) = 7rII(cc), the 
generation of topological defects in the field 9 (x) can be 
described by a term 



with Fourier transform 



P c (ey) « 2ir 2 VD 



for q -» 0. 

It is convenient to introduce the normal state tunnel 
conductance per unit of length in the SN junction^ 



Gi 



h 

2^ 



2 p w ^wP2D^2D 

2nJ ^ Z ' 



where p w (p2T>) is the normal-state local density of states 
in the wire (film), and fi w (^-2D) is the volume of the 
wire (area of the film) . Replacing the expression of the 
cooperon appearing in Eq. ([9| and noting that the resis- 
tivity in a 2D film is g 2D = [e 2 ^fr e ] _1 = [e 2 p 2 DD]~ 1 , 
we can express Eq. ^ as 



S A 



Ga 
so 



dxdx'P c (x - x') cos [9 (x) - 9 (x 1 )] , (11) 



where Ga is the dimensionless Andreev conductance in 
the SN junctiorPEH 



E 

71=1 



so 



J dx cos (2n0 (x)) , 



(12) 



where A ps = exp {— S'corc} is the "fugacity" of a phase-slip, 
and S'corc is the actio n associated with the creation of a 
single phase-slirJ 3 ° l 33 [ The term cos (2n<fi (x)) represents 
the creation of a kink of value 27rn in the 6*-field at the 
space-time point x. Assuming that X ps <C 1, we can 
neglect in the following contributions with n > 1 in S ps . 
Adding Eqs. fLT]) and ([12} we finally arrive at 



the expression of the effective phase-only action at low 
temperatures 



S — Sq + Sa + S ps , 



(13) 



describing on an equal footing the effects of fluctuation, 
dissipation and topological excitations. 



III. RESULTS 



A. Renormalization-group analysis 



G A = 



4e 2 G 2 g 2D . 



In addition to the term Sa of Eq. Q, the coupling 
t± generates contributions O (tj) and O (tjj at scales 
00 Co an d t < £o w _1 ; which renormalize the bare 



To study the properties of the model of Eq. (13 1 
at T — 0, we perform a renormalization-group analysis 
(RG) which is perturbative in the couplings Ga and A ps . 
At lowest possible order, the RG equations are found by 
performing one- and two-loop corrections in Sa and S ps , 
respectively. 



5 



We adopt a renormalization procedure that rescales 
space and time homogeneously, so as to preserve the 
Lorentz invariance of Sq. 

The renormalization of Sa involves a projection onto 
the most relevant sector of the coupling kernel P c (x), 
which is very anisotropic in space and time, obeying a 
functional RG-flow in the general case. We can simplify 
the analysis by studying different scales of interest in the 
renormalization procedure. Depending on the final scale 
A (I) ~ Lr 1 (where A (I) = Aoe 1 is the renormalized mo- 
mentum cut-off, and where Aq" 1 = £ ), we focus on the 
local part of the cooperon for A (I) < or on the non- 
local, diffusive properties for A (/) > £^ . 

We can motivate the RG analysis in the non-local 
regime by noting that the kernel P c (x) induce Josepshon 
coupling of phases over spatial distances W (r) ~ \J Dr. 
Indeed, an effective purely local kernel P° ff (x) can be ob- 
tained by integrating the spatial coordinate in Eq. (10 I, 
yielding at long times 



(x) ~ t- 2 W (t) 6 (x) for r^>D/u A 



(14) 



This approximate form is simpler to analyze, and yields 
a scaling dimension § . Note that this long-range tempo- 
ral r~ 3 / 2 differs from the standard local ohmic coupling 
r~ 2 coupling^, and further quenches fluctuations of the 
phase. A more detailed functional RG procedure involves 
an expansion of P c (x) in terms of Legendre polynomials, 
and allows to extract the scaling dimension of P c in the 
non-local limit in a systematic way (cf . Appendix [A} . 
Using that (cos (6 (x) - 9 (0))) - r' 1 / 2 *" for r -> oo (cf. 
Ref. H]), we conclude that the scaling dimension of the 
perturbative term Sa is § — 

In the local regime and for X ps — 0, our the RG- 
analysis reduces to that obtained in Ref. [T71 where de- 
tails of their derivation can be found. In this case, the 
scaling analysis of Sa is simpler to obtain, since in this 
limit P c (x) is RG-invariant, yielding a scaling dimension 
of 



We obtain the flow equations 

d J^±=G A (l)-\l s (l)K\l) 
du{l) ~ n . u{l) flW-flM 

AG a (I) 



K(l) BW 



1 



dl 



d\ ps {I) 
dl 



2K(l) 
1 



G A (l) (local), 



(15) 
(16) 

(17) 



g^^y ) G A (0 (non-local), 
= (2-K(l))X ps (l), 



where we have defined the dimensionless couplings G a = 
Gat (B' 1 ' + -B^- 1 ) and X ps = X ps VA for convenience. 
The dimensionless quantities A, B^ X,T ^> are non-universal 
and arise from the renormalization of S ps and Sa, re- 
spectively, at scales {i, ur} < A -1 (I), and are defined 
as 



A 



47T 



df f 3 e- 2KF ^F A (f) . 



A-M0 



1 1 



2 A 2 (I) '. 
1 1 



2tt 



2tt 



P c ( A" 1 (Z),</>) cos 2 ^, 



BiT) = ~2WWu? J ^^(A-^O.^-inV, 

where P c has been expressed in cylindrical coordinates 
(cf. Eq. A8l, and where f = A (I) \J x 2 + (ut) 2 and 

F\(x) — g Inf. It is interesting to point out that while 
only one parameter A arises in the rescaling of S ps due 
to space-time isotropy, the anisotropy of Sa generates 
different parameters iJ^and B^ T \ 

Note that in the local regime (P c ~ P c (r)), the prod- 
uct uK does not renormalize for X ps = 0, and thus 
B" = 0. Non-locality is thus captured by a B" > 0. 
Further, since the term Sa breaks the space-time isotropy 
within our Lorentz-invariant RG analysis (i.e., momen- 
tum shell integration homogeneous in space-time) , we ex- 
pect a renormalization of the velocity u (I), cf. Eq. ( 16 1. 
A numerical evaluation gives meaning 
that u (I) flows towards smaller values. 

As for Eq. (18 1, we note that it corresponds to the 



usual Berezinskii-Kosterlitz-Thouless (BKT) flow equa- 
tion (cf. Ref. OH] for a derivation in the context of super- 
conducting wires). 

In the limit ^Ga (0 , X ps — > 0, the properties of the 

system are dominated by the value of K (I). From Eqs. 
(17 1 and (18 1, we can define the critical values K* A = 
|) f° r t ne local (non-local) regime, and K* s = 2. 
For X ps = and K > K* A , the coupling Ga (0 flows 
towards strong coupling and eventually the perturbative 
RG analysis is no longer valid. On the other hand, for 
Ga = and K < K* s the coupling A ps (I), becomes 
relevant and eventually superconductivity is destroyed 
in the wir e, du e to the unbinding of pairs of topological 
excitation d 2 * 30 !. Note that it is not possible to determine 
the nature of the T = fixed point in this regime within 
our formalism. This issue is currently under intensive 
research^. 

Therefore, at T = and when neither Ga 

nor X ps 

vanish, the Luttingcr liquid phase is never stable, and 
the ground state is determined by a competition between 
Sa and S ps . 



B. Self-consistent Harmonic Approximation 

To further investigate the properties in the regime 
where Ga is the dominant parameter that flows to strong 
coupling, we used a self-consistent variational appr oach. 
the so-called self-consistent harmonic approximation^^. 
This method consists in finding the optimal propagator 
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9tr (q) of a harmonic (Gaussian) trial action 
2pL gtr (<?) 



that minimizes the variational free-energy 
where 

F tT = -T In / V6 e- St ' [e] . 



The minimization of the free-energy F var with respect to 
9ti (?) yields a self-consistent equation for g tr (q) 



where {e l8 ) is the value of the superconducting order 
parameter 



<e«) = 



1/4K 



1/6K 



(local) , 



(non-local) . 



(22) 



This result indicates that the order parameter develops 
long-range order, and should be compared with the case 
of isolated wires, where superconducting correlation func- 
tions follow a power-law behavior and (e' e ) — as a 
consequence of the strong quantum fluctuation^. 



C. DC transport properties 



2G a f 

9t^ (q) = 9a 1 (q) - -pt \ dx [cos - 1] 



x P c (as) exp 



— ^[l-co B (g / a;)] 5 t r (g / )J, 



(19) 



where g 1 (q) = ^w 2 ^ + ^k 2 is the propagator in the 



Luttinger liquid phase. The solutions of Eq. (19 1 read 

i , , = f 5q 1 (?) +?7Kn| (local), 
tr W I % 1 (?) + ^y^Tkl (non-local). 

(20) 

The parameter 77 is found self-consistently for the general 
case, but in the limit G a — > it reduces to 



7] = < 



2irG A £Nexp (^) 

f 2 
mi 



4AT 



8ir 2 G A VD 
T 2 



(local) , 



(non-local) . 



Note that physical solutions of the Eq. (19 1 with 77 7^ 



are found only for K > K* A . confirming the results of the 
RG analysis. In the context of the variational approach, 



it becomes clear (cf. Eq. (20 1 ) that the contribution 



of the cooperon propagator of Eq. (Ill induces ohmic 
(non-ohmic) dissipation in the local (non-local) regime. 
Evaluation of the phase-correlation function at T = 
with the optimal gtr (<z) of Eq. ( 20 1 , yields in the long 
wavelength limit 

(1 + -^L 1 (local), 
,2y/D 1 ^ 
1 H s — ^ 7-^ — (non-local), 



rjTT 2 x 2 + ADt 



(21) 



Now we address the experimentally relevant question 
of the possibility to observe some signatures of our pre- 
dictions at T = 0. To that end, we turn our attention 
to transport properties, and calculate the DC resistivity. 
We use the theoretical framework of the memory matrix, 
which is perturbative i n the p rocesses that degrade the 
current-density operator 2 * 29 * 48 ! 

, , uK2e„ n/ , 
J (as) = \76(x) . 

7T C 

Current-decay originated by phase-slips induce finite re- 
sisitivity at T < T c . At very low temperatures T«r c 
thermally activated phase slips are suppressed and re- 
sistivity is dominated by quantum phase-slips processes. 
In the absence of dissipation (Ga = 0), its expression is 
well-known and reads^ 9 -^! 



4tt 3 A 2 s Ao „ 9 / AT 



B 



- , 1 — K I cos 



(23) 



where B (x, y) is the beta function. This is a valid ex- 
pression provided that a perturbation expansion in X ps 
and Ga is possible. At finite temperatures, the effect of 
these couplings can be incorporated by replacing the bare 
parameters by the renormalized ones obtained from the 
integration of the RG-flow equations (Eqs. ( l~5fl8 )) up 
to a scakP A" 1 (I) = u {I) /2nT . 

Our results are shown in Fig. [2] where we calcu- 
late q (T) normalized to the "high-temperature" value 
To = £o/u, fixed by the short time cut-off of the the- 
ory. According with our estimations (see Sec. IV I, we 
analyze only the local regime A (/) -C ^ . We start with 
the initial conditions K (0) = 2.1,A ps (0) = 10~ 3 , for 
Ga(0) — (solid line in Fig. [2]). For comparison, we 
show the (T/To) 2K< ' ^~ 3 behavior predicted for the re- 
sistivity due to phase-s lips in t he absence of dissipation 
effects (dot-dashed line)! 29 1 30 1 33 1 . 

Interestingly, starting the RG flow with the initial val- 
ues G A (0) = 0.01, 0.05, 0.1 and 0.2, the resistivity de- 



creases faster than the (T/Tq) 



2K(Q)-3 



law corresponding 



7 



0.1 

P 0.01 
Q. 

0.001 




G A =0.00 
G A =0.01 
, G A =0.05 
G A =0.10 
G A =0.20 
(T/T ) 2K " 3 



0.001 



Figure 2: (color online) Normalized resistivity vs T/Tq. As 
the (dimensionless) Andreev conductance Ga is increased, the 
wire resistivity o(T) deviates from the law ~ y 2 ^-3 p re di c ted 
for an isolated wir e 1 1 'as a consequence of the dissipation- 
induced increase of the stiffness K (cf. Eq. (15l). 



to the isolated wire. This illustrates the stabilizing effect 
of dissipation on superconductivity, which manifests it- 
self throu gh a n increase of the stiffness K, as can be seen 
from Eq. (15 1 when parameter Ga dominates over A ps . 

Note that since the integration of the renormalization 
group flow (and consequently, the calculation of the re- 
sistivity) is perturbative, it cannot be carried beyond a 
point where either Ga (0 or X ps (I) become of order unity. 



IV. DISCUSSION AND SUMMARY 



The result of the RG flow Eqs. (|15|-(|18j) together with 
the analysis in the strong coupling regime, summarized 
in Eqs. (21 1 and (22 1, suggest that a weak coupling to 



the metallic film favors a superconducting ground-state 
with long-range order of the order parameter at T = 0, 
through a dissipation-induced quench of phase fluctua- 
tions. 

Note that this is not trivial, since a strong coupling to 
a disordered metallic film is detrimental to superconduc- 
tivity and lowers T c through the well-known proximity 
effectP. But in a low-dimensional situation at T <C T c . 
where phase fluctuations are the dominant mechanism of 
destruction of global phase-coherence, the environment 
actually favors long-range order. 

This picture is supported by experiments on disordered 
granular Pb films coated with a thin Ag metallic filnP^l, 
where it was shown that while T c decreases due to the 
proximity effect, phase stiffness actually increases at low 
enough temperatures. Also in the context of dissipative 
Josephson junctions arrays, it is well-known that the ex- 
istence of coupling to a normal metal stabilizes the su- 
perconducting phases^^. 

A similar idea has been recently suggested to produce 



an enhancement of T c in bilayered materials^, in which 
one layer has a high pairing scale but low superfluid den- 
sity while in the other layer the situation is the inverse. 
When both materials are put into contact, the T c of the 
coupled system is higher than those of the isolated layers. 
In the specific case of Luttinger liquids coupled to dissi- 
pative baths, our results are in agreement with recent 
theoretical works where the existence of superconductive 
long-range order at T = has been suggeste d 17 ! 18 !. 

In this paper, we have presented a rigorous study of a 
realistic dissipative mechanism, provided by a coupling 
to a diffusive metal, in the context of superconducting 
wires. 

Of central importance in our analysis is the cooperon 
propagator kernel P c (x — x'), which couples the field 9 
at the space-time coordinates x and x' in the dissipa- 
tive term Sa (cf. Eq. (Ill)- The physics of the ker- 
nel P c {x — x 1 ) strongly depends on the relation between 
£jv {T), the coherence length in the diffusive film, and 
the length of the wire L. Consequently, two regimes of 
interest appear: the local regime £jv (T) <C L, where 
the coupling of phases is purely local in space and non- 
local in time; and the non-local regime £jv (T) > L. 
where the phase-coupling is non-local both in space and 
time. At this point, it is interesting to determine to 
which regime actual superconducting wires would cor- 
respond. At the experimentally relevant temperature 
T ~ IK, and using typical values of D in clean metal- 
lic films^ D ~ 10 2 cm.s _1 , we obtain the estimate 
£jv (l-?0 ~ 0.1 fim. On the other hand, the temperature 
constraints to observe non-local effects can be compactly 
written as Lt (T) <C L -C £/v (T). These conditions 
require that T 3> T/vl = fi.it 2 //c^Z? (where units have 
been restored) while the length has to be kept smaller 
than £at (T). Estimating the velocity of the Mooij-Schon 
as ii ~ 10 5 m.s~ 1 for the wires in Ref. we obtain 
Tnl ~ 10 — 100 K, which exceeds the bulk-T c values es- 
timated in the range 6.9 - 7.1 K in In-Pb filmsP. The 
above estimates show that spatial non-local effects are 
elusive in actual wires (e.g., such as those studied in Ref. 
|2"T1 where L ~ 10 — 100 //m), but may eventually be ob- 
served in superconducting wires with higher T c , coupled 
to very clean substrates. 

In order t o make contact with recent transport 
experiments^^ 27 * 31 !, we have calculated the linear DC 
resistivity of a wire weakly coupled to a diffusive film, 
for different values of the Andreev conductance Ga- 
The results of Fig. [2j calculated for the local regime 
£jv (T) <C L, suggest that signatures of the predicted 
long-range order phase at T = could be observed exper- 
imentally. Indeed, since the dissipative term Sa renor- 
malizes the superconducting stiffness K to higher val- 



ues as the temperature decreases (cf. Eq. (15 1), si zable 
deviations from the predictions for an isolated wir d 29 * 30 * 
q (T) - T 2K( -°1- 3 , where K (0) is the bare stiffness, could 
be achieved at low enough temperatures. 

In summary, we have studied a thin superconducting 
wire weakly coupled to a metallic film, focusing on the 
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details of dissipation provided by the metallic cooperons 
at low temperatures. We have studied the phase diagram 
at T = within the framework of renormalization group 
and a variational harmonic approximation. We predict a 
quantum phase transition towards a superconductor with 
long-range order at T = as a function of the Andreev 
conductance G a and the superconducting stiffness K of 
the wire. Finally, we show that some signatures of this 
ordered phase could be observed in experiments of trans- 
port, manifested as an increase of the superconducting 
stiffness and consequently the exponent of g (T) ~ T v at 
low temperatures. 
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Appendix A: The cooperon propagator 

In this section we derive the expression for the 
cooperon propagator assuming weak disorder in the 
metallic film. We refer the reader to Ref. 03] for further 
details. 

When evaluating averages over the disorder potential 
inEq. (|| 

P c (r -r',T- t') = ( 52 D (r, r'; r - r') g 2 u (r, r'; r - r')) d 

a diagrammatic series (ladder diagrams) is constructed 
upon the repeated action of the Dyson's equation (in op- 
erator notation) 

92D =<72D +92D^92D 

where ' s the unperturbed electron Green's function 
in the otherwise perfect metal, and 8(r — r')V(r) = 

V r'\ is the (static) disorder potential which veri- 



r 

fies 



(V(r)) d = (Al) 
(V(r)V(r')) d = ni V 2 S(r-r') (A2) 

where rii is the concentration of impurities and V is the 
uniform component of V (r) in Fourier space. The dia- 
grammatic series in Fourier space representation is given 
by (cf. Ref. 23} 



P c (Q;n,m) 



P°(Q;n,m) 



1 - 



n,V 5 



P°(Q;n,m)' 



where 



p° (Q; n , m ) = E g2D ( p ' il>n ) g2D ~ p ' iL>n 



and where f2 is the volume of the sample. Q — k + k' rep- 
resents the center-of-mass momentum of the two electron 
system, and k and k! are the initial (i.e., before collid- 
ing with impurities) momenta of the individual electrons. 
Defining 



((Q]iv n ,iui 



ILOl = lV m — IV, 

n,V 2 



P°(Q;n,l + n) 



the real-space representation of P c reads 



Pc(r,r) 



1 



J2 ^ e -i(2v„+«,)r e i0.r x 



Q=k-\-h' n,l 



riiV 



((Q;iv n ,iu)i) 



1 - ((Q;iv n ,ioJi) ' 



(A3) 



Using the definition of the mean-free path l e = uj?r e , 
and the definition of the elastic lifetime (in Bom's 
approximation )@Sl 



— = 2 7 rp 2D n I V 2 (A4) 

we obtain in the diffusive limit (Ql e ,u>iT e <C 1) 

C {Q;iv n ,iui) »1-T e (\uji\ + DQ 2 ) , 

where we have defined the diffussion constant in d spa- 
tial dimensions D 



d 



— Replacing these re- 
sults into Eq. (|A3|, and noting that the contribution to 



£ {Q]iv n ,iuJi) is vanishingly small for v n (v n + cj;) > 0, 
we obtain 

P c (r,r)^27rp 2 ^J2 e ~ i2UnT x - 



E 



e iQ.r e -iuiiT 



DQ 2 



(A5) 



v„(v„+ui)<0 



In addition, and since the condition v n (v n + u)i) < 
must be fulfilled, we note that the limit r e — ► con- 
strains the summation over Matsubara frequencies v n to 
values near v n k, 0. Therefore, at low temperatures we 
find 



sin(|w ; |r) e 



iQ.r 



DQ 2 



We can generalize this expression to take into account 
processes that break phase coherence (magnetic fields, 
magnetic impurities, etc) 



P c (r,r) = 2p 2D 



1 

y u,,Q 



sin t) 



JQ.r 



\lji\+DQ 2 +7V 1 ' 
(A6) 
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where t v is a phenomenological phase-breaking time. At 
T = we can replace -g Y],^ 2tF Io° < ^ U3m Then, the sum 



over Q in Eq. (|A6|) yields for a 2D system 

P2D 



P c (r, t) = Re 




where Ifp ( x ) is the zeroth-order modified Besscl func- 
tion, r = \J x 1 + y 2 is the distance in the x — y plane and 
t c is an ultraviolet cutoff in time to make the integral 
in uj convergent, and which we set r c = r e . Finally we 
arrive at the expression of the cooperon in Eq. (|9| 



(A7) 



in a systematic way 



f(r,4>) = J2At (r)T e (cos < 



where the coefficients are defined as 



of 4- 1 r 1 

At (r) = — ±— J ^ d (cos 4>) f (f , 0) Pi (cos 4>) . (A12) 



Then. 



(An) 2 1 °° 

P c (r,0) = ^^-^A £ (r)^(cos0) (A13) 



Changing variables to v = cos</>, we can write Eq. (A12| 
as 



where we defined 



ie 5 « e*Df r 

P c (r, t) = Re ; 







T = T + IT P 



and where r [a, z] is the incomplete gamma function. 



Expansion in terms of Legendre polynomials 



2/ 4- 1 f 1 

Ae{r) = ^~ / duV t (u)f(f,u) 



2f+l r 1 , V £ (v) 
dv 



2 fV-i/ 2 + l' 



(A14) 



An asymptotic approximation of (f ) in the regime f — > 
oo gives 



In order to investigate the scaling dimensions of the 
bare kernel P c in the non-local regime, we can use the 
approximate expression of Eq. ( 10 1 



Pc{x) 



(AD? 



(ADr? 



and express it in terms of cylindrical coordinates 
{x = r cos 0, ut = r sin (j)} as 



P c Or) = P c (r, , 



(AD) 2 1 - 



where the definitions 

f(r,4>) 



rjj r 2 



f 2 cos 4 6 + sin 2 i 



AD Av F , 

rD = = h, 

u u 



(A8) 

(A9) 
(A10) 
(All) 



has been used and where vp is the Fermi velocity in the 
metallic film. We now expand / (f , <fi) in Legendre poly- 
nomials in order to separate radial and angular variables 



2£+l 



Jim A e (f) -^P e (0) j ^ dv V2vi + V 



2 

21+ I 



1 



Vi(0) 



'2f 



(A15) 



where from the Rodrigues formula^ we obta in Vi (0) — 
:P — ^ even (note from Eq. (|A14|) that 



Ai (f) vanish for odd values of £). Therefore, from Eqs. 
(A13land (A15) we obtain in the asymptotic limit 



Pc (r, 4>) (0) Vt (cos <f) , 



\/2 r 
7T (4L>) 2 1 



(cos</)) , 



where in the last line the expansion of the Dirac-delta 
function in terms of the Legendre polynomials has been 
used. Coming back to coordinates x and r, we obtain 



P c (x) 



nVAD 6 (x) 

V2 T§ ' 



which provides a rigourous derivation of Eq. (14 1 



10 



1 N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 

(1966) . 

2 T. Giamarchi, Quantum Physics in One Dimension (Ox- 
ford University Press, Oxford, 2004). 

3 A. O. Caldeira and A. J. Leggett, Ann. Phys. (N. Y.) 149, 
374 (1983). 

4 A. J. Bray and M. A. Moore, Phys. Rev. Lett. 49, 1545 
(1982). 

5 S. Chakravarty, Phys. Rev. Lett. 49, 681 (1982). 

6 A. Schmid, Phys. Rev. Lett. 51, 1506 (1983). 

7 G. Schon and A. D. Zaikin, Physics Re- 
ports 198, 237 (1990), ISSN 0370- 
1573, URL http : //www . sciencedirect . 
com/science/ article/B6TVP-46T4WHC-6V /2/ 
56e621d5b6c3e08daef3b06fb8b5d3c6 

8 R. Fazio and H. van der Zant, Physics Reports 355, 235 
(2001). 

9 P. A. Bobbert, R. Fazio, G. Schon, and A. D. Zaikin, Phys. 
Rev. B 45, 2294 (1992). 

10 P. Goswami and S. Chakravarty, Phys. Rev. B 73, 094516 
(2006). 

11 G. Refael, E. Demler, Y. Oreg, and D. S. Fisher, Phys. 
Rev. B 75, 014522 (2007). 

12 M. V. Feigel'man and A. I. Larkin, Chem. Phys. 235, 107 
(1998). 

13 M. V. Feigel'man, A. I. Larkin, and M. A. Skvortsov, Phys. 
Rev. B 61, 12361 (2000). 

14 B. Spivak, A. Zyuzin, and M. Hruska, Phys. Rev. B 64, 
132502 (2001). 

15 M. V. Feigel'man, A. I. Larkin, and M. A. Skvortsov, Phys. 
Rev. Lett. 86, 1869 (2001). 

16 A. H. Castro Neto, C. de C. Chamon, and C. Nayak, Phys. 
Rev. Lett. 79, 4629 (1997). 

17 M. A. Cazalilla, F. Sols, and F. Guinea, Phys. Rev. Lett. 
97, 076401 (2006). 

18 S. N. Artemenko and T. Nattermann, Phys. Rev. Lett. 99, 
256401 (2007). 

19 W. Little, Phys. Rev. 156, 396 (1967). 

20 J. S. Langer and V. Ambegaokar, Phys. Rev. 164, 498 

(1967) . 

21 D. E. McCumber and B. I. Halperin, Phys. Rev. B 1, 1054 
(1970). 

22 N. Giordano, Phys. Rev. Lett. 61, 2137 (1988). 

23 N. Giordano, Physica B 203, 460 (1994). 

24 C. N. Lau, N. Markovic, M. Bockrath, A. Bezryadin, and 
M. Tinkham, Phys. Rev. Lett. 87, 217003 (2001). 

25 M. Tian, J. Wang, J. S. Kurtz, Y. Liu, M. H. Chan, T. S. 
Mayer, andT. E. Mallouk, Phys. Rev. B 71, 104521 (2005). 

26 M. Zgirski, K.-P. Riikonen, V. Touboltsev, and K. Aru- 
tyunov, Nano Lett. 5, 1029 (2005). 

27 F. Altomare, A. M. Chang, M. R. Melloch, Y. Hong, and 



C. W. Tu, Phys. Rev. Lett. 97, 017001 (2006). 

28 A. T. Bollinger, R. C. Dinsmore III, A. Rogachev, and 
A. Bezryadin, Phys. Rev. Lett. 101, 227003 (2008). 

29 T. Giamarchi, Phys. Rev. B 46, 342 (1992). 

30 A. D. Zaikin, D. S. Golubev, A. van Otterlo, and G. T. 
Zimanyi, Phys. Rev. Lett. 78, 1552 (1997). 

31 A. Bezryadin, C. N. Lau, and M. Tinkham, Nature 404, 
971 (2000). 

32 H. P. Biichler, V. B. Geshkenbein, and G. Blatter, Phys. 
Rev. Lett. 92, 067007 (2004). 

33 K. Y. Arutyunov, D. S. Golubev, and A. D. Zaikin, Phys. 
Rep. 464, 1 (2008). 

34 J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 

(1973) . 

35 F. W. Hekking and Y. V. Nazarov, Phys. Rev. Lett. 71, 
1625 (1993). 

36 F. W. Hekking and Y. V. Nazarov, Phys. Rev. B 49, 6847 
(1994). 

37 I. J. R. Aitchison, P. Ao, D. J. Thouless, and X.-M. Zhu, 
Phys. Rev. B 51, 6531 (1995). 

38 A. van Otterlo, D. S. Golubev, A. Zaikin, and G. Blatter, 
Eur. Phys. J. B 10, 131 (1999). 

39 A. L. Fetter and J. D. Walecka, Quantum theory of many- 
particle systems (McGraw-Hill, New York, 1971). 

40 K. B. Efetov and A. I. Larkin, Sov. Phys. JETP 39, 1129 

(1974) . 

41 J. E. Mooij and G. Schon, Phys. Rev. Lett. 55, 114 (1985). 

42 C. Bruder, R. Fazio, and G. Schon, Physica B 203, 240 
(1994). 

43 E. Akkermans and G. Montambaux, 



le me- 

soscopique des electrons et des photons (EDP Sci- 
ences/CNRS, Paris, 2004). 

44 M. Tinkham, Introduction to Superconductivity (McGraw- 
Hill, 1996), 2nd ed. 

45 A. M. Lobos, A. Iucci, M. Miiller, and T. Giamarchi, (in 
preparation). 

46 P. G. De Gennes, Superconductivity of Metals and Alloys 
(W. A. Benjamin, New York, 1966). 

47 R. P. Feynman, Statistical Mechanics (Benjamin, Reading, 
MA, 1972). 

48 W. Gotze and P. Wolfle, Phys. Rev. B 6, 1226 (1972). 

49 L. Merchant, J. Ostrick, R. P. Barber, and R. C. Dynes, 
Phys. Rev. B 63, 134508 (2001). 

50 E. Berg, D. Orgad, and S. A. Kivelson, Phys. Rev. B 78, 
094509 (2008). 

51 A. Morpurgo, private communication. 

52 M. Abramowitz and I. A. Stegun, Handbook of mathemat- 
ical functions : with formulas, graphs and mathematical 
tables (Dover, New York, 1965). 



